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We study metric-afRne gravity (MAG) inspired cosmological models. Those models were statis- 
tically estimated using the SNIa data. We also use the cosmic microwave background observations 
and the big-bang nucleosynthesis analysis to constrain the density parameter f2^,o which is related 
to the non-Riemannian structure of the underlying spacetime. We argue that while the models 
are statistically admissible from the SNIa analysis, complementary stricter limits obtained from the 
CMB and BBN indicate that the models with density parameters with a a 6 scaling behaviour 
are virtually ruled out. If we assume the validity of the particular MAG based cosmological model 
throughout all stages of the Universe, the parameter estimates from the CMB and BBN put a 
stronger limit, in comparison to the SNIa data, on the presence of non-Riemannian structures at 
low redshifts. 
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Tj" ■ I. INTRODUCTION 

(N . 

Astronomical observations brought important changes in modern cosmology While the type la supernovae 
(SNIa) data are most often employed, the cosmic microwave background (CMB) observations and big-bang nucle- 
osynthesis (BBN) analysis can also be used. They allow the exotic physics of cosmological models to be checked 
against the observational data [2| • There is an increasing effort to develop some cosmological and astrophysical tools 
to search for new physics beyond the standard model. 

The metric-afRne gravity (MAG) cosmological model with the Robertson- Walker symmetry was investigated by 
three different groups. First, it was considered the model with triplet ansatz in vacuum Second, it was considered 
the dilational hyperfluid model £|. Third, it was shown that on the level of the field equations the special case of the 
MAG model is equivalent to a model in the Weyl-Cartan spacetime if we choose a model parameter in the special 
form (a6 = —04) ; 5j. Moreover after redefinition of some variables the second and third approach gives the same set 
of dynamical equations. The analysis of constraints on parameters in the MAG model can be addressed in all three 
approach, but we adopt the last one proposed by Puetzfeld and Chen 

Note that in the model with dust matter on the brane, apart from dark radiation which scales like a -4 , there is a 
correction of the type a -6 to the Einstein equations on the brane which arise from the influence of a bulk geometr 

The term scaling like a~ 6 also appears in the Friedmann- Robertson- Walker model with spinning fluid 
It is possible to establish formally the one to one correspondence between the MAG model and either the Randall- 
Sundrum brane model when positive values of the non-Riemannian contribution to effective energy is admitted or the 
spinning fluid filled cosmology when this contribution is negative. However, if one takes the pure Randall-Sundrum 
type model then there is a constraint on the brane tension parameter coming from the theory itself. The brane tension 
parameter is not less than about (100 GeV) 4 . The MAG model is free from such a theoretical constraint. 

In our further discussion we examine the flat models which is motivated by the CMB WMAP observations an d 
consider the following formula for the Friedmann first integral 
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TABLE I: Best fit estimation of the model parameters from the SNIa data 



priors 




r2 r ,o "*/j,o 




M 


9 

X 


"V,o ^ 0; i2 m ,o > 





0.14 0.012 


0.848 


15.945 


175.75 


^V,o > 0; fi m ,o = 0.3 




0.005 


0.695 


15.965 


177.30 


fim,0 > 





0.14 0.012 


0.848 


15.945 


175.75 


^m,0 = 0.3 




0.005 


0.695 


15.965 


177.30 


Qm,o > 0; Q r ,o = 0.0001 


0.16 


- 0.029 


0.811 


15.945 


175.97 


n m , = 0.3; fi r , = 0.0001 




- 0.005 


0.695 


15.965 


177.30 



TABLE II: Maximum likelihood estimation of the model parameters with 2a errors from the SNIa data 



priors 



^V,o > 0; Qm,o > 
Qif>,o > 0; fl m ,o = 0.3 
n m , > 

"m,0 ~ 0.3 

fim.o > 0; f2 r , = 0.0001 
fi m ,o = 0.3; Q r ,o = 0.0001 



m,0 



fit 



Oa.o 



g+0.24 






q+0.15 00 q+o.04i 
0.00 u -0.00 u - uuy -0.009 
q+0.017 



^0.04 
-0.00 



-0.000 



q+0.27 +0.21 onf) q+0.042 
u -0.00 u -0.00 u - uuy - 







,+0.09 



0.14 H 



-0.036 

o.oo2±S:gg 

0.0281°;™ 
0.005+°°}? 



o.82o±g:i°g 

0.680+°,;™ 



0.800 



+0.110 



0.6801°;°*° 



0.810 
0.695 



+0.130 
0.150 
+0.015 
0.O15 



where H — dlna/dt is the Hubble function, t is the cosmological time, a = a(t) is the scale factor, ft m .o, ^A,o 
and fi^o are the density parameters for dust matter, the cosmological constant and fictitious fluid which mimics 
"non-Ricmannian effects" , respectively. Their values in the present epoch are marked by the index "0" . All density 
parameters satisfy the constraint condition 

^m,0 + ^A,0 + ^r,0 + ty),0 = 1- (2) 

The density parameter for the fictitious fluid is defined as p| 



where 



k 2 f 1 3a 



144a V h 

The sign of the parameter v is undetermined and it can assume both positive and negative values. 



II. CONSTRAINT FROM THE SNIA 



Let us start from the reestimation of the models parameters by using the latest sample of SNIa data [llj. The 
motivation to study the SN constraint is to find the best estimation available from the latest data which gives the 
narrowest constraint for this method. In the next sections we compare it with constraints obtained from other 
methods. 

Riess et al.'s sample contains 157 type la supernovae [llj. We consider the flat model with and without priors on 
the VL^fi and Jim, o- We assume that the former can be of any value or only nonnegative, and the latter is nonnegative 
or equal 0.3 [12]. We estimate the best fits of the model parameters (Table QJ. Additionally we find the maximum 
likelihood estimates with the errors at 2a level (Table [p. 

We find that that the estimates of the parameter Sl$ t o are very close to zero although positive apart of one case 
when it is zero. We can conclude that the estimate of this parameter is order of magnitude of 0.01. 

To illustrate the results of the maximum likelihood analysis of the model we draw the levels of confidence on 
Figure Q 

The MAG model fits well to SNIa data. We consider the model with any value of fl rj o we obtain the value of £l m .o 
equal zero as best fit and maximum likelihood estimator, while fixing the small amount radiation (f2 ri o |f 3| ) gives the 
low density matter universe. The estimation of the Hubble constant gives the value close to 65 km/s MPc. 
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FIG. 1: The contours with lcr and 2cr confidence levels for fi,/,,o versus Ho, £Ia,o, fim.o, an d ^r,o from the SNIa data. 



III. CMB PEAKS IN THE MAG MODEL 



The hotter and colder spots in the CMB can interpreted as acoustic oscillation in the primeval plasma during the 
last scattering. Peaks in the power spectrum correspond to maximum density of the wave. In the Legendrc multipole 
space these peaks correspond to the angle subtended by the sound horizon at the last scattering. Further peaks answer 
to higher harmonics of the principal oscillations. 

It is very interesting that locations of these peaks are very sensitive to the variations in the model parameters. 
Therefore, it can be used as another way to constrain the parameters of cosmological models. 

The acoustic scale Ia which puts the locations of the peaks is defined as 



£ A = n 



Jo 



Zdcc dz' 



H(z') 



r 

J z, 



dz' 



where 



h{z) = H y/n mfi (i + z)3 + n r . (i + zY + n^(i + zf + n Afi 

and c s is the speed of sound in the plasma given by 



dpcs 



[n 7i0 (l + z) + 6^, (l + z) 3 



s dp cS 314,0 + 40^(1 + z) + 6fl^ (l + z) 3 ' 
Knowing the acoustic scale we can determine the location of m-th peak 



(4) 



(5) 



(6) 



£m ~ lA(m - 4> m ) 



(7) 
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TABLE III: Values of Q^ t o an d location of first three peaks 



Hubble constant 










H = 65 km/s MPc 


3 x 10" 11 


220 


522 


825 




7 x 1(T 14 


220 


523 


826 




-1.4 x 10" 10 


223 


530 


847 


H = 72 km/s MPc 


3.7 x 10" 11 


220 


522 


823 







220 


521 


821 




-1.3 x 10" 10 


224 


530 


847 



where (j) m is the phase shift caused by the plasma driving effect. Assuming that f2 m ,o = 0.3, on the surface of last 
scattering z^ cc it is given by 



0.267 



K^dec) 



0.3 



0.267 



1 PrfZdcc) 



0.3 p m (>dec) 



0.267 



1 ^ ri0 (l + z doc ) 



0.3 



0.3 



o.i 



(8) 



where f2b,oft 2 = 0.02, r(zd cc ) = Pr^decVPm^dcc) = ^r,o(l + Zdcc)/^m,o is the ratio of radiation to matter densities at 
the surface of last scattering. 

The CMB temperature angular power spectrum provides the locations of the first two peaks [Til and the 
BOOMERanG measurements give the location of the third peak They values with uncertainties on the level la 
are the following 



h = 220.1^;*, t 2 = 546tig, 4 = 845±^. 



Using the WMAP data only, Spergel et al. [14( obtained that the Hubble constant Hq = 72 km/s MPc (or the 
parameter ft = 0.72), the baryonic matter density f2b,o = 0.024ft -2 , and the matter density r2 tllj o = 0.14/i~ 2 which 
give a good agreement with the observation of position of the first peak. 

To find whether cosmological models give these observational locations of peaks we fix some model parameters. Let 
the baryonic matter density f2b,o — 0.05, the spectral index for initial density perturbations n — 1, and the radiation 
density parameter 



fir.o = ^7,0 + £Vo = 2.48/i~ 2 x 10~ 5 



1.7ft" 2 x 10~ 5 = 4.18ft" 2 x 1CT 5 



0) 



which is a sum of the photon f2 7i o and neutrino f2„.o densities. 

Assuming rZ m o = 0.3 and ft = 0.72 we obtain for the standard ACDM cosmological model the following positions 
of peaks 



?i = 220, l 2 = 521, 



= 821 



with the phase shift <fi m given by JSJ • 

From the SNIa data analysis it was found that the Hubble constant has lower value. Assuming that Hq = 65 km/s 
MPc (or ft = 0.65), we have f2 rj o = 9.89 x 10 -5 from equation ©. In further calculation we take Sl rj o = 0.0001. If we 
consider the standard ACDM model, with fi m ^ = 0.3, fib.o = 0.05, the spectral index for initial density perturbations 
n=l, and h = 0.65, where sound can propagate in baryonic matter and photons we obtain the following locations of 
first three peaks 



£1 = 225, i 2 = 535, 



= 847. 



We find some discrepancy between the observational and theoretical results in this case. Now it is interesting to check 
whether the presence of the fictitious fluid fl^o change the locations of the peaks. 

The properties of the fictitious fluid f^o are unknown. In particular, we do not know whether the sound can or 
cannot propagate in this fluid. But we assume that sound can propagate in it as well as in baryonic matter and 
photons. We consider both values of the Hubble constant and assume that h = 0.65 or ft = 0.72. The results of 
calculations of peak locations and the values of the parameter fi^o are presented in Table ITTT1 

If we choose the Hq = 65 km/s MPc then we obtain the agreement with the observation of the location of the 
first peak for three non-zero values of the parameter O^,.o- As it is shown on Figure [5] there are two positive and one 
negative values of this parameter for which the MAG model is admissible. 

All these distinguished values of f^o are m agreement with the result obtained from SNIa because the 2a confidence 
interval for this parameter obtained from the SNIa data contains these three points. While the SNIa estimations give 
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FIG. 2: The location of the first peak in function of f2^,,o- 



the possibility that 0^,0 is equal zero, the CMB calculations seem to exclude this case because the zero value of fi^o 
requires the first peak location at 225. 

If we choose the Ho = 72 km/s MPc than one of positive values of f^o move to zero, while the second one move a 
little to the right. 

We also calculated the age of the Universe in the MAG model. We find that the difference in the age of the Universe 
is smaller than I mln years for all three values of Sl^o- Assuming that f2 m; o = 0.3 the age of the Universe is 14.496 
Gyr for Hq = 65 km/s MPc, and 13.088 Gyr for Hq = 72 km/s MPc. The globular cluster analysis indicated that the 
age of the Universe is 13.4 Gyr |l7|. 

IV. CONSTRAINT FROM THE BBN 

It is well known that the big-bang nucleosynthesis (BBN) is the very well tested area of cosmology and does not 
allow for any significant deviation from the standard expansion law apart from very early times (i.e., before the onset 
of BBN). The prediction of standard BBN is in well agreement with observations of abundance of light elements. 
Therefore, all nonstandard terms added to the Friedmann equation should give only negligible small modifications 
during the BBN epoch to have the nucleosynthesis process unchanged. 

In our opinion the consistency with BBN is a crucial issue in the MAG models where the nonstandard term a~ 6 in 
the Friedmann equation is added (see also discussion in 0])- This additional term scales like (1 + z) 6 . It is clear that 
such a term has either accelerated (f^o > 0) or decelerated (fii/,,o < 0) impact on the Universe expansion. Going 
backwards in time this term would become dominant at some redshift. If it would happen before the BBN epoch, the 
radiation domination would never occur and the all BBN predictions would be lost. 

The domination of the fictitious fluid fl^ should end before the BBN epoch starts otherwise the nucleosynthesis 
process would be dramatically modified. If we assume that the BBN result are preserved in the MAG models we 
obtain another constraint on the amount of f^o- Let us assume that the model modification is negligible small during 
the BBN epoch and the nucleosynthesis process is unchanged. It means that the contribution of the MAG term Sl^o 
cannot dominate over the radiation term f2 rj o ~ 1CU 4 before the beginning of BBN (z ~ 10 s ) 

0^,0(1 + z) 6 < n 1%0 {i + z) 4 =► |^, |<io- 20 . 

The values of fl^, o oc 1CP 2 obtained as best fits in the SNIa data analysis as well as the smallest nonzero value 
of fl^o — 7 x 1CU 14 calculated in the CMB analysis arc unrealistic in the light of the above result. If we take into 
consideration the maximum likelihood analysis of SNIa data we have the possibility that the value of Q^o is lower 
than 1 1 20 1 in the 2a confidence interval. In the case of the CMB analysis only the value of the Hubble constant 
close to 72 km/s MPc gives the zero or close to zero value of fi^Q. 
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V. CONCLUSION 



The paper discusses observational constraint on "energy contributions" arising in certain cosmological models based 
on MAG. In particular it is focused on the nonstandard term a~ 6 . We test this model against the SNIa data, the 
location of the peaks of the CMB power spectrum, and constraints from the BBN. 

The MAG model fits well to SNIa data and the estimations give the amount of fluid fl^o to be order of magnitude 
0.01, and the Hubble constant is close to 65 km/s MPc. Let us note that these results are compatible with constraints 
from FRIIb radio galaxies and X-ray gas mass fractions 0] . 

The CMB analysis gives that the Hubble constant is 72 km/s MPc which gives the too low age of the universe 
in comparison with the age of globular clusters. Taking lower value of the Hubble constant obtained from SNIa 
estimation resolves the problem of the age. However, the location of the first peak shifts to the right and is in conflict 
with the observed location. The introducing of the non-Riemannian structure of the underlying spacetime moves the 
location of the first peak back and this MAG model agrees with the CMB observations. The analysis of the CMB in 
this model cannot distinguish the character of the fictitious fluid and we do not know whether the parameter fi^o is 
positive or negative. 

The absolute values of f2^.o obtained in the MAG model from the CMB analysis with h — 0.65 seems to be too 
large in comparison to the limit obtained from the BBN analysis. Using the BBN analysis we pointed out that the 
MAG part of the energy density to its present density parameter is of order 10~ 20 . The limit of this order leads to 
the conclusion that the MAG model is virtually ruled out. However, we must remember that we insist that the MAG 
model does not change the physics during and after the BBN epoch. In this context, the merit of the SNIa analysis 
is its independency from any assumption on physical processes in the early Universe. 
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